home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / slaed5.z / slaed5
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. SSSSLLLLAAAAEEEEDDDD5555((((3333FFFF))))                                                          SSSSLLLLAAAAEEEEDDDD5555((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      SLAED5 - subroutine computes the I-th eigenvalue of a symmetric rank-one
  10.      modification of a 2-by-2 diagonal matrix   diag( D ) + RHO  The diagonal
  11.      elements in the array D are assumed to satisfy   D(i) < D(j) for i < j
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      SUBROUTINE SLAED5( I, D, Z, DELTA, RHO, DLAM )
  15.  
  16.          INTEGER        I
  17.  
  18.          REAL           DLAM, RHO
  19.  
  20.          REAL           D( 2 ), DELTA( 2 ), Z( 2 )
  21.  
  22. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  23.      This subroutine computes the I-th eigenvalue of a symmetric rank-one
  24.      modification of a 2-by-2 diagonal matrix
  25.  
  26.      We also assume RHO > 0 and that the Euclidean norm of the vector Z is
  27.      one.
  28.  
  29.  
  30. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  31.      I      (input) INTEGER
  32.             The index of the eigenvalue to be computed.  I = 1 or I = 2.
  33.  
  34.      D      (input) REAL array, dimension (2)
  35.             The original eigenvalues.  We assume D(1) < D(2).
  36.  
  37.      Z      (input) REAL array, dimension (2)
  38.             The components of the updating vector.
  39.  
  40.      DELTA  (output) REAL array, dimension (2)
  41.             The vector DELTA contains the information necessary to construct
  42.             the eigenvectors.
  43.  
  44.      RHO    (input) REAL
  45.             The scalar in the symmetric updating formula.
  46.  
  47.      DLAM   (output) REAL
  48.             The computed lambda_I, the I-th updated eigenvalue.
  49.  
  50.  
  51.  
  52.  
  53.  
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.